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Abstract 

This article derives the first two moments of the two versions of the Riesz distribu- 
tion in the terms of their characteristic functions. 



1 Introduction 

There is no doubt about the important role played by Wishart distribution in the con- 
text of multivariate stati stics and random matrix theory. Based in the Riesz measure, 
Hassairi and Laimil (| 2 lh proposed a generalisation of the Wishart distribution, which they 
termed Ri esz distribution. Some of their main properties and rela ted d istributions have bee n 
studied b ylHassairi and Laimil (|200ll ). lFaraut and Koranvil (|1994l ) and lHassairi et all ([20051 ) . 
Recently, iDiaz-Garcial (|2012n proposes two versions of the Riesz distribution for real normed 
division algebras and some of their properties are also being studied. 

In particular, the characteristic function was obtained for both versions of the Riesz 
distribution, but a topic that has been disregarded, is the study of its moments. As indicated 
in the conclusions section, in particular, these moments can be used to study the asymptotic 
normality of the Riesz distribution. 

This article studies the first two moments for the two versions of the Riesz distribu- 
tion. Section [2] reviews some definitions and notations on the matrix algebra and some 
special functions with matrix arguments on real symmetric cones, also, are summarised the 
two definitions of the Riesz distribution, their corresponding characteristic functions and 
the section is finalised with obtaining the first two matrix derivatives of the characteristic 
functions. The main results are proposed in Section |3l 



2 Preliminary results 

2.1 Matrix algebra and special function with matrix argument 

A detail ed discussion of theory of matri ces a nd special f unctio n with matrix argument can be 
found in lMagnus and Neudeckeil ()l988f ) and lMuirheadl (|l982[ ). respectively. For convenience, 
we shall introduce some notation, although in general, we adhere to standard notation forms. 

Let 2l mxn be the set of all m x n matrices over SR. Let A e 2l mxn , then A' G 2t" xm 
denotes the transpose. It is denoted by & m the real vector space of all S G 2t" ixm sucn that 
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S = S . In addition, let be the cone of positive definite matrices S £ 2l mx «\ Thus, <p m 
consist of all matrices S = X X, with X G 2l™ xm ; then *p m is an open subset of & m . 

r m [a] denotes the multivariate Gamma function for the space G m . This can be obtained 
as a particular case of the generalised gamma function of weight k for the space 6 m with 
ft = (fci, k2, ■ ■ . , k m ), ki > k2 > ■ ■ ■ > k m > 0, taking ft = (0, 0, . . . , 0) and which for 
Re(a) > (m — l)/2 — k m is defined by (see Gross and Richards! . Il987 ). 



r m [a,ft] = / etr{-A}|A|°-( Tn - 1 V 2 - 1 9K (A)(dA) 



i(m — 



W[[r[a + h-(i-l)/2} 



= [a] K T m [a], 

where etr(-) = exp(tr(-)), | • | denotes the determinant, and for Ae6, 



q K (A) = |A m | fe ™ J] |A, 



ki-kt 



(1) 



(2) 



(3) 



with A„ = (a rs ), r, s = 1, 2, . . . ,p, p = 1, 2, . . . , m is termed the highest weight vector, see 
Gross and Richards! (|l987t ). Also, 



T m [a] = [ etrl-AHAr^- 1 )/ 2 - 1 ^) 

_ n ^-^Y[T[a- (i-l)/2], 
i=i 

and Re(a) > (m - l)/2. 

In other branches of mathematics the highest weight vector q K (A.) is also terme d the 
generalised power o f A and is denoted as A K (A), see iFaraut and Koranvi (1994) and 



Hassairi and Laimil (|2001l ). 



Additional properties of q K (A), which are immediate consequences of the definition of 
q K (A) and the following property 1, are: 



1. if Ai, . . . , A m , are the eigenvalues of A, then: 

m 



(4) 



g K (A- 1 )= (Z - 1 (A)= (Z _ K (A), 

3. if k = (p, . . . ,p), then: 

<Zk(A) - |A|* 
in particular if p — 0, then q K (A) = 1. 

4. if t = {tx,t 2 ,...,t m ), ti >t 2 > ■■■ > t m > 0, then: 

q K+T (A) = q K (A)q T (A), 
in particular if r = (p,p, ■ ■ ■ ,p), then: 

q K+T {A) = q K+p {A) = \A\ p q K {A). 



(5) 
(6) 

(7) 

(8) 
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5. Finally, for B € 21™*™ such that C = B*B € e m , 

q K (BAB*) = q K (C)q K (A) 

and 

^(B^AB*- 1 ) = (< ZK (C))- 1 g K (A). 



(9) 



(10) 



Remark 2.1. Let P(6 m ) denote the algebra of all polynomial functions on & m , and 
'Pki&m) the subspace of homogeneous polynomials of degree k and let V K '(& m ) be an irre- 
ducible subspace of V{& m ) such that: 

V k (& m ) =^0^(6™)- 

K 

No te that q K is a homogeneous polynomial of degree k, moreover q K G P K (& m ), for references 
see 



Gross and Richards! (Il987l) 



In @, [a] K denotes the generalised Pochhammer symbol of weight k, defined as: 

m 

wk = n( a -( i_i )/ 2 ) fc * 

i=l 

m 

^mim-V/i JJ P [ a + fc . _ ( ? _ X )/ 2 ] 



r m [a] 



r m [a, k] 



r m [a] 

where Re(a) > (m — l)/2 — fe m and 

(a)! = a(a + 1) ■ • • (a + i — 1), 

is the standard Pochhammer symbol. 

An altern atively definition of the generalised gamma function of weight k is proposed by 
Khatril (Il966h . which is defined as: 



r m [o,-/s] = / etr{-A}|A|°-( m - 1 )/ 2 - 1 g K (A- 1 )(dA) 

•/Apse™ 



_ ro(m-l)/4 



(11) 



J|r[a-*s - (m-i)/2] 



(-l) fe r m [a] 



[-a+ (m- l)/2 + l]„' 



(12) 



where Re(a) > (m — l)/ 2 + fci . 



Also recall that from iMagnus and Neudecken (|l988f ): 
1. Let A be an m x n matrix and B a p x q matrix. The mp x nq matrix is defined as: 

aiiB • • ■ ai„B 

^mlB • • • <2 mn B 

and is termed the Kronecker product of A and B and written A ® B. 
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2. Let A be an m x n matrix and Aj its j-th column, the vec A is the mn x 1 vector 



vec A 



A 2 



3. tr AB = vec' A' vec B. 

4. tr ABCD = vec' D'(C <g> A) vccB = vec' D(A' ® C) vecB'. 

5. For any matrix X = (xij) £ 5ft mx ™, dX denotes the matrix of differentials, dX 

(tfajy). 

6. In particulafl if F : U mxm -> 5R mxm of rank m, then: 

(a) d|F(X)|* = |F(X)|nr(^(X))- 1 dF(X), 

(b) dF(X)" 1 = -F(X)- 1 dF(X) J F(X)- 1 . Also 

(c) dAXB = AdXB, and 

(d) d tr AXB = tr AdXB. 

7. If / : K mx " -> vec'X = (vecX)' and dvec'X = (dvecX)' then: 

(a) if dvec/(X) = vec'BdvecX, its obtained that: 

<9vec/(X) ,„ , df(X) „ 

» V Y = vec'B and -^i = B. 
a vec' X aX 

(b) if d 2 vec/ (X) = dvec'X B dvecX then: 

9 2 vec/(X) _ 1 



tV = o(b + b'). 



d vec X<9 vec'X 2 
Finally, note that for a matrix A of order m x m, their submatrices 

an • • • air, 



Ai = an, A 2 = 



an a i2 
a 2 i a 2 2 



■ * ■ A m — 



a?7ii 



can be written as: 



A,; = E'AE,; 



(13) 



where the columns of Ej : m x i are the first i columns of the identity matrix I m and 
E< = [Ij ]' is such that 

i X m — z 

(a) E^E 4 = I 4 

(b) E+ = Ej, hence, E+ = E<, where A + denotes the Moore- Penrose inverse of A, 

(c) E m Im: 

(d) and 





m — ixi 



ft . o J 





i x m — i 

o 



m — ixi m — ixm — i 



1 These results exist for more general conditions, sec Magnus and Ncudcckcr (1988, Theorem 1, pp. 149) 
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Closely related to the differentiation are the commutation matrix K mn and the matrix 

The commutation matrix K m „ is such that for a to x n matrix A, 

K mn vec A = vec A' . 

with, K^ ln = = K nm . If to = rt, is written as K m instead of K mn . The main property 
of this matrix is: Let A be an m x n matrix, B a p x q matrix. Then: 

K pm (A <g> B) = (B ® A)K gn . 

By other hand, the matrix N m is such that: 

N — N' — N 2 — K N — NK 

Note that for A an m x m symmetric matrix, 

vec A = N m vec A. (14) 

Finally consider the following matrix factorisation. 

Proposition 2.1. If A. is a non-negative definite m x m matrix then there exist a non- 
negative definite to x to matrix, written as A 1 / 2 such that A = A 1 / 2 A 1 / 2 . 

2.2 Riesz distributions 

This section summarise the densities and the ir correspond i ng ch aracteristic functions for the 
two versions of the Riesz distribution. From iDiaz-Garcfal (|2012h . 

Definition 2.1. Let £ E *p m and k = (kx, fe, • . . , k m ), ki > k% > ■ ■ ■ > k m > 0. 

1. Then it is said that X has a Riesz distribution of type I if its density function is: 

^^^'^""•W" ,15) 

for X G *}3 m and Rc(a) > (m — l)/2 — fc m ; denoting this fact as X ~ £H^(a, /c, S). 

2. Then it is said that X has a Riesz distribution of type II if its density function is: 

r r^L,, etr{-^X}\X\ a -( m -^ 2 -%(X-')(dX) (16) 

for X e *p m and Re(a) > (to — l)/2+ kx', denoting this fact as X ~ (a, k, £). 
From Diaz-Garcial (|2012h and using (J7]) it is obtained that: 
Lemma 2.1. Let £ s 93 m a^rf k = &2> • • • j ^m), fci > &2 > • • • > k m > 0. 

1. TTierc if X ~ 9i^(a, k, S) ?is characteristic function is: 

X (T) = g K+a ((l m -lE^TE 1 /*)"^ (17) 
/or Re(a) > (m - l)/2 - fc m . 

TTien i/ X ~ (a, «, X) iis characteristic function is: 

X (T) = g K _ a (l m - zX^TX 1 / 2 ) (18) 
/or Re(a) > (to - l)/2 + fc x . 
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2.3 Differentiation 

Finally, consider the following result about differentiation. 

Lemma 2.2. Let £ G ty m and k — (k±, fe, • • • , k m ), k\ > > ■ • ■ > k m > 0. Then if 
X ~ £H^j(a, k, S) it is obtained: 

1. 

dvec X (T) 



dvec'T 



I?, and 




yfot - ti+i) vec' A 

i=l 

<9 2 vec (/> x (T) „ [A, . 9 . . . 

^ xV / = i 2 N m { V (U -ti+i 2 vecA,vec'A I 

' vec Ta vec' T ' 

k 2=1 

m 

( <l - vecAi vec' A 

X! X! ( <l _ - ij+i) vec Aj vec' A 



f 

+ 2 



jr(U - *<+i) (A, ® A,) In to<Zt ((l m - iS 1 / 2 TS 1 / 2 ) 1 j 



where 



A a = S 1 / 2 E Q (E^(l m -zS 1 / 2 TS 1 / 2 )E Q ) 'e^S 1 / 2 , a = i,j, 
r = ac + a = (fi, . . . ,t m ) and Re(a) > (m — l)/2 — fc m . 

Proof. By 6(a), 6(b) and Lemma 12.11 defining k + a = t = (t\, . . . ,t m ), with t m +i = 0, 
E m = I m and 



A Q = S 1 / 2 E Q (e; (l m - i S 1 / 2 TS 1 / 2 ) E Q ) 1 E^S 1 / 2 , a 



i,3, 



it is obtained that 

<# x (T)=i 
Hence, vectorising, 

d vec </> x (T) = i 



^(U-t i+1 )tTAidT 



jZ(u ~ ^+0 vec ' a » 



T ( (I m - iE^TE 1 / 2 ) ) dvecT 



from where applying 7(a) and (|14[) . the desired result is obtained. 
Analogously, differentiating again, 

{m 
Y2(U - i l+ i) 2 tr A, ; dTtr AidT 
i=i 

m 

+ 5^(tj - ti+xjtr A^TA^T L T (l m - zS 1 / 2 TS 1/2 
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And vectorising using 4, it is got; 



d 2 vec0 x (T) = i 2 i - t i+ i) 2 dvec'TvecAivec' AidvecT 



. i=l 



+ y~^y~] (U —U+i)(tj — tj + x)dvec' TvecAi vec' AjdvecT 



- t l+1 )dvec' T (A, ® A,) dvecT \q T ( (l m - iS 1/2 TS 1/2 



From where, by applying 7(b) and (fT4")l , the outcome sought is obtained. □ □ 

Analogously, one has: 

Lemma 2.3. Let £ 6 <p m and k = (fci, k2, . . . , k m ), k\ > k.2 > • • • > k m > 0. Then if 
X ~ (a, k, X) it is obtained: 



d vec X (T) 
dvec'T 



N m ? T fl m - iS 1 / 2 TS 1/2 ) , 



,9. and 



<9 2 vec0 x (T) 
SvecTdvec'T 



y^fti - U+i) vec' Aj 



N m < y^fe - < l+ i) 2 vec Aj vec' A* 

m 

X] X (** - ) ~ vec Ai vec ' A 

1 

/] (*» ~ - vecAj vec' A 

- ti+i) (A, <g> A 4 ) lN m g T (l m - zE 1 / 2 TS 1 / 2 



where: 



A a = ^ 2 E a (e; (l m - zS^^TS 1 / 2 ) E a ) 1 K* 1/2 
t = K — a = (ii, . . . , t m ) a^rf Re(a) > (m — l)/2 + k±. 

Proof. This is a verbatim copy of the proof of Lemma 12.21 □ 



a = i, J, 



3 Moments of Riesz distributions 

This section proposed the main result. 

Theorem 3.1. Let X S ty m and n — (fci, fe, ■ ■ • , fc m ), fci > &2 > 
1. Then if X /ias a iiiesz distribution of type I, 



■>k m >0. 



□ 
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(a) E(X) = (k m + o)E + ( k * - Zsi-OS^EiEjEVa. 

i=l 

(7>^ and Cov(vecX) is 

7/7 1 

(fc m + a)N m (S ®-£)+J2(ki- ki-i)N m (s^E^i/a 8 £Va E . E ^. s i/a 

»=i 

/or Re(a) > (m - l)/2 - fe ro . 

Then if X /ias a Riesz distribution of type II, 

m— 1 

(a) E(X) = -(km - a)S - - fc^iJS 1 /^,^^ 1 / 2 . 

(7>^ and Cov(vecX) is 

m—l 

-{k m - a)N m (S ® £) _ (fc, - fci_i)N m (sV2 EiE / s i/2 g, sVa^EjS 1 /: 

8=1 

/or Re(a) > (m- l)/2 + k 1 . 
Proof. Results are immediately from lemmas 12.21 and 12.31 remembering that: 

dvec (/> x (T) 



E(vec X) 



i <9vecX 



T=0 



and 
where 



Cov(vecX) = E(vecXvec'X) - E(vccX) E(vec' X) 

d 2 vec X (T) 



(19) 



E(vecXvec'X) = 



T=0 



i 2 d vec X<9 vec' X 
In order to ensure that Cov(vecX) = Cov'(vecX) it is necessary that 

E(vecX) E'(vecX) 

be a symmetric matrix. Then, proceeding as in the case of E(vecXvec' X) = (B + B')/2, 
consider the following equivalent definition; 

Cov(vecX) = E(vecXvec'X) - i{E(vecX) E(vec' X) + [E(vecX) E(vec'X)]'}, 

which alternative definition, coincides with (|19[) when E(vec X) E'(vec X) is a symmetric 
matrix and at the same time ensuring that Cov(vec X) is a symmetric matrix. □ □ 

Observe that in Theorem 13.11 1 and 13.11 2. are defined as a — n/2, X — > 2S and k = 
(0, . . . , 0) the Wishart case is obtained. Moreover, 

1. E(X) =nS, 

2. Cov(vecX) = 2nN ro (S ® S) = n{l rn 2 + K m )(S ® £), 
see lMuirheadl (|l982l p. 90) and lMagnus and Neudeckeil (|l988l p. 253). 
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4 Conclusions 



Observe that if Yi, ■ • • , Yjv, are independent p-dimensional random vectors, with N > p, 
such that the random matrix X, is defined as: 



A' 



X = ^Y. i Y.; = Y'Y, Y' = (Y 1 ,---,Y JV ) 



has a Riesz distribution type I, then by multivariate central limit theorem. iMuirheadl (|1982 , 
p. 15) and Theorem 13. 11 if 

1 N 

N 1 



and 



1 N 

S(n) = -V(Y i -Y)(Y i -Y) / , TV = n + 1 



the asymptotic distribution as n — > oc of 



,1/2 



vec S(n) vec 



£+ V (fc ' ^ vecS^^E^S 1 / 2 



is: 



Af m * ( ,^±__N m (S®E)+^ 



m — 1 



N m f 53V2E E ' S i/2 g S i/2 EiE ' S i/2 



Note that this asymptotic multivariate normal distribution is singular, moreover, its rank 
is m(m + l)/2. Also, observe th a t if a = n/2, S — > 2S and k = (0, ... ,0) the asymptotic 
result is obtained by iMuirhead ( 1982 . pp. 90-91) for the Wishart case. The author is 
currently studying in detail the distribution of the random matrix Y and some of their basic 
properties. 
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